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WEAK INCIDENT SHOCK INTERACTIONS WITH 


MACH 8 LAMINAR BOUNDARY LAYERS 

By Louis G. Kaufman n 1 and Charles B. Johnson 
Langley Research Center 

SUMMARY 

Weak shock-wave interactions with boundary layers on a flat plate were investigated 
experimentally in the Langley Mach 8 variable-density tunnel for plate-length Reynolds 
numbers ranging from 0.46 x 10® to 2.5 x 10®. The undisturbed boundary layers were 
laminar over the entire plate length. Pressure and heat- transfer distributions were 
obtained for wedge-generated incident shock waves that resulted in pressure rises rang- 
ing from 1.36 to 4.46 (both nonseparated and separated boundary- layer flows). The 
resulting heat-transfer amplifications ranged from 1.45 to 14. The distributions fol- 
lowed established trends for nonseparated flows, for incipient separation, and for lami- 
nar free-interaction pressure rises. The experimental results corroborated established 
trends for the extent of the pressure rise and for certain peak heat-transfer correlations. 
Because of the many factors that strongly affect laminar boundary layers, and because 
transition frequently occurred prior to the end of the interaction region, the data did not 
support the simpler correlations of peak heating with peak pressure that have been pro- 
posed previously. 


INTRODUCTION 

Shock waves that impinge on a surface can greatly amplify the local heat transfer 
and pressure loads on the surface. The adverse pressure gradient associated with a 
sufficiently strong shockwave causes a laminar boundary layer to separate from the 
surface. Although the problem is well recognized and has received much attention 
(refs. 1 to 4), there are still no adequate simple methods for predicting pertinent geo- 
metric features of the separation region and the increased heating at reattachment. 

In general, when boundary -layer separation occurs because of shock interaction, 
there is a region of reverse flow (fig. 1) that terminates where the separated boundary 
layer reattaches to the surface, with attendant high pressures and heating rates. The 
compression waves emanating from the separation location coalesce into a shock, and 

1 Research Department, Grumman Aerospace Corporation; work performed as a 
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the pressure rises to a "plateau" value. For laminar boundary layers, the heating rate 
drops below the undisturbed value at the onset of separation. Immediately downstream 
of the compression waves are expansion waves that emanate from the top of the separa- 
tion bubble. Downstream of the expansion waves there are recompression waves (coa- 
lescing into a shock) that occur with reattachment of the separated flow, and the pressures 
and heating rates can attain very large values. Therefore, obtaining meaningful experi- 
mental pressure and heat-transfer data in a shock— laminar-boundary-layer interaction 
region involves careful design of experimental tests and model. For example, in the 
present tests to minimize extraneous effects, a simple two-dimensional configuration 
was chosen: wedge -generated shock waves interacting with boundary layers on a sharp- 
leading-edge flat plate. The wedges were designed to provide as long a region (stream- 
wise) of constant pressure as possible in the given wind tunnel. As sketched in figure 1, 
the wedges were positioned as far forward as possible. The wedge trailing edge was 
positioned near the plate bow wave to minimize reflected waves. The experiments were 
conducted in the Langley Mach 8 variable -density tunnel for plate length Reynolds numbers 
ranging from 0.46 x 10® to 2.5 x 10®. 

The experiments provided pressure and heating- rate distributions for incident shock 
interactions with initially laminar boundary layers, and schlieren photographs of the inter- 
action flow field. Many investigators (refs. 5 to 11) examined the peak heating at reattach- 
ment and several proposed simple correlations in the form h pk °c (peak heating- 

proportional to the peak pressure raised to a power n). However, the values of the expo- 
nent n in these current correlations vary from 0.5 to 1.3, an unacceptable difference. 

An objective of the subject work was to reexamine these correlations in an attempt to 
resolve their discrepancies. In addition to this objective, the present data are compared 
with existing, simple, analytical methods for predicting "free interaction" pressure rises 
(refs. 12 and 13), plateau pressure levels (refs. 2 and 14 to 16), criteria for incipient 
separation (refs. 17 to 21), and the extent of laminar separation (refs. 21 to 25). 

SYMBOLS 

c specific heat of model material, j/kg-K 

Cf coefficient of friction of boundary layer 

Cp specific heat of air flow, J/kg-K 

F(s) similarity function for free-interaction pressure rise (see eq. (6)) 

H ratio of disturbed to undisturbed heat-transfer coefficients at same location 
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S = fX 


a 

y 


heat- transfer coefficient, w/ m^K 
length of pressure rise upstream of x^ 
Mach number 


Prandtl number 


Stanton number, shy^pCpVj 


exponent in correlations of peak heating with peak pressure 

(^/^)disturbed 


ratio of disturbed to undisturbed pressure ratios, 


pressure, N/m^ 


( p / p t) u „ dl 


undisturbed 


heating rate, W/m2 

Reynolds number based on free- stream conditions and a unit length of 1 meter 


Reynolds number based on local undisturbed flow conditions and length Xj 


recovery factor 
- *o) 

temperature, K 
time, sec 
velocity, m/sec 

streamwise distance measured from flat-plate leading edge, m 


reference location used for free-interaction pressure rise similarity 
function F, x = x^ at F = 0.81 


turning angle for incipient 


separation, radians 


ratio of specific heats (taken as 1.4 herein) 
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Ap uncertainty in measured pressure, N/m2 

At' change in Prandti-Meyer angle for isentropic compression, deg 

6 boundary -Layer thickness, m 

Q angle of attack of wedge shock generator, deg 

A proportionality constant in equation (8), 

I u viscosity, N-sec/m2 

p density, kg/m 3 

r thin-wall thickness, m 

flow-deflection angle at reattachment 

Subscripts : 

aw adiabatic wall conditions 

e conditions at outer edge of boundary layer 

i incident shock location on plate surface in absence of boundary layer 

incip conditions for incipient separation of laminar boundary layer 

m model material (stainless steel) 

plat plateau value 

pk peak (maximum) value 

sep local conditions at separation location 

si shear layer 

t total (stagnation) conditions 
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w wall conditions 

0 local undisturbed flow conditions at location where interaction pressure 

starts to rise 

1 local undisturbed flow conditions 

2 conditions immediately downstream of incident shock 

3 conditions downstream of reflected incident shock 

00 free -stream conditions 

A comma is used to indicate a double subscript. For example, Hp pk refers to 
the value of H at P pk . 


APPARATUS AND TEST PROCEDURES 
Model Design 

The model configurations and tunnel flow conditions were chosen to provide data 
for shock -wave boundary -layer interactions with a minimum of extraneous effects. A 
simple two-dimensional configuration was chosen: wedge -generated shock waves imping- 
ing on a flat -plate surface. The wedge and plate were designed to keep extraneous waves 
(such as those from the wedge trailing edge) as far away from the interaction region as 
possible. (See fig. 1.) Small wedge angles were used to insure obtaining small inter- 
action flow regions and primarily to insure obtaining fully reattached flow ahead of the 
influence of any extraneous effects. Tunnel flow conditions were chosen to yield laminar - 
boundary -layer flow over the entire flat -plate surface in the absence of an incident shock 
wave. 

Photographs of the model (fig. 2) show the wedge shock generator sting mounted to 
a vertical strut. The instrumented flat -plate shock receiver is cantilevered forward from 
the same strut. The wedge has a span of 15.25 cm; the flat plate has a span of 17.46 cm. 
The chords of the wedge and flat plate are indicated in the profile view sketched in fig- 
ure 3. The instrumented center line of the flat plate is well within the region bounded by 
the tip Mach cones. 

The wedge shock generator could be set at various streamwise stations and at vari- 
ous heights above the flat -plate shock receiver. In order to avoid reflections of the plate 
bow wave from the wedge surface, the wedge trailing edge is positioned just above the 
calculated position of the bow wave from the plate leading edge. (See fig. 1.) Indeed, it 
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was hoped that this bow wave would partially cancel the expansion fan emanating from the 
wedge trailing edge. 

The wedge shock generator can be set at any one of three incidence angles with 
respect to the free-stream flow: 6 = 1°, 3°, or 5°. The effective wedge angles are some- 

what different for two reasons. The displacement thickness growth of the boundary layer 
on the wedge surface increases the effective wedge angle, whereas the interaction with the 
plate bow wave decreases the effective wedge angle. The maximum static -pressure rise 
on the plate surface is that associated with the shock caused by the effective wedge angle 
as well as its reflection from the flat -plate surface. If these considerations are taken 
into account, the chosen angles are sufficient to provide data for unseparated flows as 
well as for small regions of separated flow. (See refs. 19 to 21.) Larger angles would 
cause extensive regions of separated flow and increase the probability of the free shear 
layer becoming transitional prior to reattachment (ref. 26) and/or the reattachment region 
being influenced by the wave from the wedge trailing edge. Thus, the shock generator is 
designed to yield as large a region of ''clean” shock impingement interaction flow as 
possible. 

There are two (interchangeable) flat -plate shock receivers. Both have machined 
sharp leading edges. One is instrumented with 48 pressure orifices evenly spaced along 
the plate center line; and one, with thermocouples at the same locations. (See fig. 3.) 

The first instrumentation location is 13.94 cm downstream of the leading edge; the pres- 
sure orifices and thermocouples are evenly spaced 0.508 cm apart. The iron-constantan 
thermocouples are spot welded to the inner surface of a "thin skin” section of the plate 
used for heat-transfer measurements. The skin thickness in this region is 0.076 cm. 

The pressure instrumented plate has a constant thickness of 0.318 cm. 

Test Conditions 

The experiments were conducted in the Langley Mach 8 variable density tunnel. 

This is a blowdown wind tunnel with a closed circular test section, with a diameter of 
45.72 cm, that can provide Mach 8 free-stream flows for unit Reynolds numbers R 
varying from 0.60 x 10® to 30 x 10® per meter. (See refs. 27 and 28.) The tunnel has 
a model injection system directly beneath the test section. Although full injection to 
tunnel center line requires approximately 0.8 second, the model passes through the tun- 
nel wall shear layer in just 0.2 second. Therefore, the model experiences nonuniform 
flow heating for just 0.2 second. A 12. 5-m -diameter sphere and an 18.3-m-diameter 
sphere were used in conjunction to provide a vacuum system with sufficient capacity to 
enable the low pressure measurements to stabilize within the tunnel running time (approx- 
imately 4 minutes). 
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The experiments were conducted at four tunnel stagnation pressure p t levels. 
They are listed in table I along with the values of the total temperatures T^.. The cali- 
brated free-stream Mach numbers corresponding to these {relatively low) tunnel 

pressure levels are listed {ref. 27) as well as the values of the free-stream unit Reynolds 
numbers per meter. Pressure, heat transfer, and schlieren flow photographic data were 
obtained for all model configurations (6 = 1°, 3°, and 5°). 

TABLE I. - TUNNEL FLOW CONDITIONS 


Pt 

T t ,K 


R/m 

MN/m 2 

psia 

0.48 

70 

710 

7.68 

1.3 x 106 

.71 

103 

710 

7.73 ! 

1.9 

1.51 

219 

730 

7.80 | 

3.7 

2.96 

429 

1 

770 

7.87 

6.5 


Test Procedures and Data Reduction 

For the pressure and heat -transfer experiments, the tunnel flow was first started and 
established, the model was injected into the tunnel flow, data were recorded, and then the 
model was retracted prior to stopping the tunnel flow. Pressures were measured by using 

39 multirange capacitance -type transducers and 9 single-range strain-gage transducers. 
Heat-transfer data were obtained by using thermocouples and the thin-wall transient- 
temperature technique. Schlieren flow photographs were taken during each pressure 
and heat -transfer tunnel run. Additional tunnel runs were made to obtain schlieren flow 
photographs with the model mounted further forward in the test section. The model wall 
temperature T w was measured and recorded after each pressure run was completed. 

Although 80 percent of the pressure orifices are connected to multirange 
capacitance -type transducers, which are much more accurate than the single-range 
strain-gage transducers, the accuracy of the pressure gages is considered herein to be 
that of the strain-gage transducers; that is, Ap = ±8.6 N/m 2 (0.00125 psi). The result- 
ing uncertainties in terms of free-stream static pressures vary from ±0. lSp^ for the 
lowest pressure runs to iO.OSp^ for the highest pressure runs. 2 

The temperature time history of each thermocouple is recorded at the rate of 

40 readings/second on magnetic tape by an analog -digital data-recording system. These 
time histories are used to obtain the temperature rise rates (dT W/ /dt), and these rates are 

2 This is the calibrated accuracy of the pressure instrumentation. The scatter in 
the pressure data, as evidenced in the following section, is considerably larger. 
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used to obtain the heating rates q from 

dT w 

a -Oct — — 

4 w 'm m m dt 


( 1 ) 


where o and c are the density and specific heat of the model material (stainless 
m ™ 

steel), and r m is the thickness of the thin wall (0.076 cm). 

Once the model starts to experience heating, the error due to conduction of heat 
along the model surface increases with time, whereas the aerodynamic heating rate 
decreases with time. Thus, the percent of heating due to conduction relative to the indi- 
cated aerodynamic heating increases rapidly with time. Therefore, the data must be 
reduced as early as possible after aerodynamic heating starts. 

For these experiments, the temperature -time derivative in equation (1) is obtained 
by least-squares fitting a linear temperature-time variation to the temperatures recorded 
during the interval from 0.1 to 1.1 seconds after the model reaches the tunnel center line. 
The simple linear variation is best suited for reducing the scatter in the relatively low 
heating-rate data. The scatter results from the electronic "noise" in the data recording 
system. The short time intervals and the relatively low heating rates lead to small vari- 
ations in wall temperature along the model surface, and negligible conduction effects. 

The experimental heat-transfer coefficient h is calculated from 


h = 



(2) 


where 


T aw = T . + r ( T t ' T 4 < 5 > 

The laminar boundary-layer recovery factor ^r ~ ^Np r ~ 0.85j is used to calculate the 
adiabatic wall temperature T aw for the reduction of all the heating data herein. 

The heat-transfer data are obtained during the initial heating of the model while the 
wall is relatively cool (T w « 0.4T^. The longer run times required for the pressure 
readings to stabilize result in the pressures being recorded while the model is somewhat 
hotter (T w ~ 0.5T|J. Although this difference is small in terms of the total temperature, 
it is undesirable because wall temperature affects both boundary- layer transition and 
separation. (See ref. 14.) This fact should be kept in mind when examining and compar- 
ing the pressure and heating rate distributions. 

As noted previously, the model is injected into the tunnel flow for the pressure and 
heat-transfer runs. The tunnel windows do not extend the full length of the test section 
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immediately above the injection cabin. When the injection system is used, only the for- 
ward section of the model is visible through the tunnel windows. In order to obtain 
schlieren flow photographs of the interaction region, the model has to be mounted directly 
in the forward part of the tunnel test section, 18 cm upstream of the injection system port. 
The model wall temperatures for these schlieren photograph runs are similar to those for 
the pressure runs (T w ~ 0.5T t ). In order to show the entire interaction flow region, it is 
necessary to combine schlieren photographs from two different tunnel runs, one being the 
tunnel run to obtain the plate pressures and the other being the tunnel run to get schlieren 
photographs of the interaction region with the model mounted directly in the tunnel. 

EXPERIMENTAL RESULTS 


Undisturbed Flows 

In order to obtain "undisturbed” pressure and heating rate distributions on the flat- 
plate surface, the wedge shock generator and sting are disconnected from the strut (see 
fig. 2) and removed. The flat plates are then tested alone to obtain the pressure and heat- 
ing rate distributions undisturbed by any incident shock waves . These distributions are 
used in determining undisturbed reference conditions for the pressure rises and heating- 
rate amplifications caused by the generated shock-wave— boundary-layer interactions. 

Curves faired through the measured undisturbed pressure distributions, nondimen- 
sionalized with respect to the total pressure for each tunnel run, are plotted in figure 4. 
The pressures on the plate surface p^ are 20 to 60 percent larger than the free-stream 
static pressures p . These increased pressures are caused by the growth of the 
boundary-layer displacement thickness; they agree within the data scatter shown in fig- 
ure 4 with other Mach 8 experimental data (ref. 9) and also with weak viscous interaction 
theory (ref. 29, p. 349). Unfortunately, the undisturbed pressures did not stabilize during 
the tunnel run for the lowest total pressure ^p = 480 000 N/m2j, and therefore are not 
plotted. Corresponding to the increased pressures on the plate surface, the local flow 
Mach numbers M^ are less than the free-stream Mach numbers listed in table I. 

The undisturbed heating rate distributions, expressed in terms of Stanton number 
Ngj., are plotted in figure 5. Herein, the Stanton number is 


N St = 



(4) 


where Cp, p, and V are the specific heat, density, and velocity of the local flow at the 
outer edge of the boundary layer. The distributions are characteristic of laminar- 
boundary-layer heating distributions on flat plates and are compared with theoretical 
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distributions calculated by using an implicit finite-difference scheme of Anderson and 
Lewis (ref. 30) for laminar boundary layers for the faired experimental pressure distri- 
butions. (See fig. 5.) 


Interaction Flows 

Salient aspects of the pressure and heat-transfer data for the interaction flows are 
listed in table II for all wedge shock generator angles. The incident shock locations 
obtained from schlieren flow photographs, are the same for the pressure and heat-transfer 
runs. The unit free-stream Reynolds numbers per meter are listed for each pressure 
and heat-transfer tunnel run with the shock generator attached. For the pressure runs 
the total pressure, the peak pressure rise ratio P pk = P p k/ P i’ the locat;ion the P eak 
pressure x pk , and the length f of the pressure rise upstream of x i are indicated. 

The upstream extent of the pressure rise is nondimensionalized with respect to the cal- 
culated (ref. 30) boundary- layer thickness for the undisturbed flow 6 Q at the location 
where the interaction flow pressure rise starts. (The distance between pressure taps is 
comparable to the boundary- layer thickness calculated at these locations.) The last two 
columns list the undisturbed local flow Mach number at the location x 0 and the Reynolds 
number based on the distance x Q and local flow conditions. For the heat-transfer runs, 
the total temperature, the ratio of model wall temperature to total temperature, the peak 
heating-rate amplification H pk s h pk /h lf the location of the peak heating x pk , and the 
value of the peak heating rate h pk are indicated. 

Schlieren flow photographs for the 1° wedge interactions are shown in figure 6. As 
described in the last paragraph of the preceding section, it was necessary to combine 
schlieren photographs from two different tunnel runs in order to show the entire interac- 
tion flow region. Detailed observation of the photographs reveals that the wedge -gene rated 
shocks initially are stronger than inviscid flow shocks for a 1° wedge, but are bent down- 
stream (weakened) by the bow waves from the plate leading edge. Extending the part of 
the wedge shock downstream of the plate bow wave to where it would strike the plate sur- 
face in the absence of a boundary layer locates what is referred to herein as the incident 
shock location on the plate surface, x^. (See fig. 6, E/m = 1.1 x 106 .) 

The pressures measured on the plate surface for these interactions are nondimen- 
sionalized with respect to the free-stream total pressure for each tunnel run. The inter- 
action pressure distribution is then taken as the ratio of p/ for the disturbed flow to 
p/Pj. for the undisturbed (flat- plate) flow: 


P = 


( p / p t ), 

( p / p t), 


disturbed 


undisturbed 


( 5 ) 
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TABLE II.- SALIENT EXPERIMENTAL RESULTS FOR INTERACTION FLOW 


(a) Pressure data 


*• 

deg 

X P 

m 

R/m 

p t’ 

MN/m 2 

P pk 

x pk» 

111 

l /\ 

M o 

R x,o 

1 

0.253 

1.4 x 106 

0.53 

— 

— 

-- 

— 



.260 

2.0 

.74 

1.57 

0.332 

7 

7.03 

0.442 x 10 6 


.269 

3.7 

1.52 

1.55 

.343 

11 

7.39 

.818 


.278 

6.3 

2.86 

1.50 

.302 

14 

7.62 

1.487 

3 

0.253 

1.3 x 10® 

0.48 

— 

— 


— 



.264 

2.0 

.77 

2.10 

0.317 

10 

7.03 

0.410 x 10 6 


.275 

4.0 

1.63 

2.17 

.322 

15 

7.39 

.844 


.283 

6.6 

3.00 

2.35 

.332 

22 

7.62 

1.459 

5 

0.266 

1.4 x 10 6 

0.50 

— 

— 


— 



.271 

1.9 

.71 

3.30 

0.327 

13 

7.03 

0.342 X 10 6 

I 

.278 

3.7 

1.50 

3.84 

.327 

20 

7.39 

.722 


.283 

6.4 

2.90 

4.46 

.353 

28 

7.62 

1.280 


(b) Heat-transfer data 



x i> 

m 

R/m 

T t> 

K 

T w/ T t 

H P k 


— 

8 


0.253 

1.3 X 106 

1 

0.412 

1.67 

0.271 

4.67 


.260 

1.8 

1 1 

.417 

1.51 

.292 

4.46 


.269 

3.6 

| | 

.396 

1.58 

.312 

6.33 


.278 

6.6 

775 

.379 

1.45 

.312 

9.05 

3 


1.2 x 10 6 

Wm 

0.414 


0.302 

6.96 


.264 

1.8 

El 

.415 


.322 

7.18 


.275 

3.9 

737 

.400 


.327 

10.42 


.283 

6.4 

767 

.385 


.368 

15.89 

5 

0.266 

1.3 x 10 6 

730 

0.405 

4.55 

0.302 

11.55 


.271 

2.0 

724 

.407 

6.75 

.368 

12.91 


.278 

3.6 

748 

.396 

9.98 

.373 

24.85 


.283 

6.5 

766 

.386 

14.05 

.373 

45.32 











This ratio minimizes effects caused by minor variations in p t from run to run. The 
undisturbed values of pyp^ are taken from the curves faired through the data (fig. 4) in 
order to avoid multiplying the data scatter. 

Interaction pressure distributions for the 1° shock generator are plotted in figure 7. 
The unit Reynolds numbers shown are those for the tunnel runs with the shock generator 
installed. (The unit Reynolds numbers for the corresponding undisturbed flow tunnel runs 
are shown in fig. 4.) For the smallest Reynolds number run (R/m = 1.4 x 10 6 , see 
table II), the pressures had not stabilized and are not plotted. The plotted pressure 
ratios start to rise (at x Q ) considerably upstream of the incident shock locations. The 
incident shock locations, obtained from examination of the schlieren flow photographs, 
are indicated by the vertical lines with the appropriate symbol denoting the correspond- 
ing unit Reynolds number. The length of the pressure rise upstream of Xj is referred 
to as "f” herein. 


Calculated values of the overall pressure rise ratios (P 3 /P 1 ) for inviscid flows are 
indicated by the horizontal dashed lines in figure 7. These values are obtained in the 
following manner. The measured undisturbed-flow pressures Pj/P^. are used to deter- 
mine the local Mach number Mj of the flow over the flat-plate surface in region 1. 

(See sketch of inviscid flow included in fig. 7.) The incident shock wave angle (down- 
stream of the plate bow wave) is measured from the schlieren flow photographs. This 
angle and Mj define the effective deflection of the local flow (from oblique shock rela- 
tions). Oblique shock relations are then employed in calculating the pressure rise across 
the incident shock pg^/pj and the local Mach number Mg in region 2. If an inviscid 
reflection of the incident shock from the plate surface (equal flow deflection angle) is 


assumed, the pressure rise across the reflected shock P^Pg calculated for Mg 
from oblique shock relations. The pressure ratios Pg/P^ anc * P 3 /P 2 are mu ^P^ ec ^ 
to obtain the overall pressure rise Pg/P^ across the incident shock and its reflection 
from the plate surface in the absence of boundary- layer effects. 


The peak pressure rise ratios listed in table II are chosen at the locations where 
the pressure distributions appear to reach their maximum values, exclusive of trailing - 
edge effects. For example, in figure 7, the three furthest downstream pressure ratio 
values are excluded. Examination of the corresponding schlieren photographs (fig. 6 ) 
reveals the possibility of a compression wave reflecting from the wedge boundary layer 
and influencing the pressures in the extreme downstream part of the instrumented section 
of the flat plate. 


Both the undisturbed (flat plate) and disturbed flow (interaction) heat-transfer coef- 
ficient distributions are plotted in figure 8 for the 1° shock generator. For the lowest 
unit Reynolds number flow, the heat-transfer coefficient distribution remains essentially 
unchanged by the interaction until just upstream of the incident shock location. It then 
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rises to a peak value downstream of the incident shock location. Holden (ref. 31) notes 
that this behavior is characteristic of fully attached, low Reynolds number, relatively 
thick boundary -layer flows with pressure rises. For the higher Reynolds number flows, 
the interaction heating rates drop below the undisturbed flow values upstream of the 
incident shock location. This decrease is associated with the thickening of the laminar 
boundary layer caused by the interaction. The heating increases to a peak value down- 
stream of the incident shock location. The initial decrease and then sharp increase 
in the interaction heating distribution for the highest Reynolds number flow is represent- 
ative of laminar boundary- layer flows approaching separation (ref. 31). 

Schlieren flow photographs of the interaction flows for the 3° wedge shock genera- 
tors are shown in figure 9. The unit Reynolds numbers shown in the schlieren photo- 
graphs are those corresponding to the photographic tunnel runs; they vary by less than 
10 percent from the corresponding unit Reynolds numbers for the pressure and heat- 
transfer tunnel runs. (See table II.) 

Interaction pressure ratio distributions for the 3° wedge shock generators are 
plotted in figure 10. Similar to the distributions for 6 = 1°, the pressures start to rise 
well upstream of the incident shock locations. The distributions do not exhibit any pla- 
teau regions of nearly constant pressures, but rather show a continuous rise to their 
maximum values, which is indicative of attached laminar boundary- layer flows approach- 
ing separation. (See ref. 31.) The measured pressure rises closely approximate the 
calculated values of Pg/Pj- 

The undisturbed and interaction flow' heat-transfer rate coefficient distributions for 
0 = 3° are presented in figure 11. There is a very small region of decreased heating 
for the lowest Reynolds number interaction. The flow can be classified as fully attached. 
(See ref. 31.) Particularly evident in the distribution for the highest Reynolds number, 
the heat-transfer rates continue to increase dowmstream of the interaction region. This 
effect can be attributed to the shock interaction disturbing the boundary layer and causing 
transition. 

Data for the largest angle wedge shock generators (0 = 5°) are presented in fig- 
ures 12 to 14. The schlieren flow' photographs (fig. 12) were scaled to obtain the incident 
shock angles and the incident shock locations indicated in figure 13. All the pressure dis- 
tributions rise to values that are fairly well approximated by the calculated values of 

*V p i' 

The interaction heat-transfer coefficient distributions (see fig. 14) are more 
rounded in the region of decreased heating than for the 0 = 3° distributions, which 
exhibit a more cusped shape. Holden (ref. 31) notes that cusped shapes are indicative 
of flows approaching incipient separation, whereas separated flows exhibit rounded dis- 
tributions of reduced heating rates. The steady downstream rise of the heating rate 
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distributions indicates the onset of boundary -layer transition prompted by the disturbance 
caused by the incident shock wave (ref. 31). Although the Reynolds numbers based on the 
total instrumented plate length (x = 0.38 m) are less than 2.5 x 10 6 , the disturbance caused 
by an incident shock wave enhances the onset of boundary -layer transition, particularly 
for free shear layers. (See ref. 26). 

THEORETICAL ANALYSES, COMPARISONS, AND CORRELATIONS 


Free Interaction Pressure Rises 

Although the maximum pressure rise and the upstream extent of the pressure rise 
depend strongly on the shock strength, the initial pressure rise should be independent of 
the mechanism generating the shock wave. This concept of a ’’free interaction” (defined 
by Chapman, Kuehn, and Larson (ref. 32)) between the viscid boundary layer and the invis- 
cid external flow' in the upstream portion of the interaction region has been well estab- 
lished (refs. 12 and 13). 

The pressure ratio distributions from where the ratios start to rise above P = 1 
to where the ratios begin to diverge for the different wedge shock generator angles are 
plotted in figure 15. The pressure ratios are plotted against (x - x o yx 0 , where the 
subscript o indicates undisturbed conditions at the upstream location where the dis- 
turbance pressure rise begins. The initial pressure ratio rises are shifted so that they 
all start at x = x Q , and therefore initially follow the same general trend. Of course, the 
stronger shock generator data approach larger pressure ratio values than do the weaker 
shock generator data. No consistent Reynolds number trends are apparent. 

The effects of the different shock generator angles are eliminated partially by plot- 
ting the data in terms of the universal similarity function F(s) suggested by Carriere, 
Sirieix, and Solignac (ref. 13): 


F(s) = 


2(P - 1) Ay 
yM 0 2 C f,o 


( 6 ) 


where s s (x - x Q ) 7(x f - x Q ), Ay is the change in the Prandtl- Meyer angle for isentropic 
compression Ay = f(P), y is the ratio of specific heats (y = 1.4 herein), c { is the 
coefficient of friction, subscript o indicates undisturbed conditions at x = x Q (the loca- 
tion where the pressure rise starts), and x^ is a reference location. The location Xj, 
which scales the extent of the pressure rise, is determined by requiring that F = 0.81 at 
x = Xj. This empirical value (F = 0.81 at s = 1.00) has been found (ref. 13) to correlate 
separated laminar boundary -layer pressure rises well for a wide range of test conditions. 
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In figure 15, the pressure ratio distributions are forced through P*1 at x = x 0 ; 
in figure 16, the distributions are forced through F = 0.81 at x = x^. Although the sim- 
ilarity function F does reduce effects of different shock generator angles, there remains 
an inconsistent Reynolds number effect. 

Plateau and Incipient Separation Pressure Rises 

As sketched in figure 1, pressure distributions corresponding to shock-induced 
laminar-boundary-layer separation frequently display a "plateau” region wherein the 
pressure varies only slightly. Many investigators have postulated that the increase in 
pressure to the plateau value should be independent of the mechanism causing separation, 
and that the plateau pressure coefficient should scale with Reynolds number. (See refs. 2 
and 14 to 16.) Extensive results from many experiments have been correlated (refs. 2 
and 15) to obtain the following expression for the laminar plateau pressure rise: 



In the preceding expression, R sep is the Reynolds number based on distance from the 
leading edge to the location of separation. For the current experiments, R sep is based 
on local flow conditions and the distance x Q . Although there is appreciable scatter in the 
experimental plateau pressure levels (standard deviation ±10 percent about eq. (7)), the 
expression is representative of laminar plateau pressure levels for a wide range of test 
conditions . 

Equation (7) is used to calculate the plateau pressure rises shown in table III for 
the present experiments. Within the accuracy of equation (7) (±10 percent), the calculated 

TABLE IE.- CALCULATED PLATEAU PRESSURE RATIOS 



plateau pressure rises are the same for all three shock generators {6 = 1°, 3°, and 5°) for 
each nominal value of the unit Reynolds number. (The plateau pressure is a "weak" func- 
tion of R sep , which is taken as R x o herein. The values of the fourth root of R xo 
agree within 10 percent for 9 - 1°, 3°, and 5° for each nominal value of the unit 
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Reynolds number; see table III). In using equation (7), the location of the pressure rise 
R sep must be known before P plat can be calculated. 

If a laminar boundary layer is well separated, then there will be a region of nearly 
constant pressure (plateau pressure) that can be estimated by using equation (7). The 
minimum pressure rise required to cause laminar boundary- layer separation (incipient 
separation) is referred to as the incipient separation pressure rise (Pi nC i p )- Values of 
p incip can be calculated b y using oblique shock relations for Mj and an incipient 
deflection angle C4 nc jp given by (refs. 17 to 20): 


MjO 


incip 


XMj 3 


( 8 ) 


The proportionality constant X, as given by Ball (ref. 17), is approximately unity for 
T w * 0.5T t . As noted above, T w ~ 0.5T t for the pressure data; therefore, X = 1 is 
used in evaluating a inci p (radians) from equation (8). Calculated values of the pres- 
sure rise ratios required to cause incipient separation are listed in table IV for the four 
nominal values of the unit Reynolds numbers for the present experiments. The measured 


TABLE IV.- CALCULATED PRESSURE RATIOS 
FOR INCIPIENT SEPARATION 


R/m 

P. - 
incip 

1.3 x 10 6 

2.2 

1.9 

1.9 

3.7 

1.7 

6.5 

1.6 


peak pressure rises for the 0 = 3° and 5° interaction flows (see table II) exceed the cal- 
culated values of P inc ^ p . However, in those present cases where separation does occur, 
the extent of the reverse flow region is too small to produce a substantial region of nearly 
constant "plateau” pressures at the calculated levels indicated in table III. (See figs. 10 
and 13.) 


Upstream Extent of Pressure Rise 

The upstream extent of a laminar pressure rise depends on several parameters: 
the overall strength of the pressure rise, the boundary-layer thickness, wall temperature, 
the character of the boundary layer at reattachment, and others. (See refs. 21 to 25.) 
Although there are no adequate empirical methods for reliably predicting the upstream 
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extent of laminar pressure rises for a wide range of conditions (ref. 2), the effects of 
various parameters on laminar pressure rises can be described at least qualitatively. 

(See refs. 21 to 25.) 

The undisturbed boundary-layer thickness at the start of the pressure rise 6 Q 
intuitively appears to be appropriate as a scaling length for the upstream extent of the 
pressure rise caused by an incident shock wave. Values of £ /6 Q from the present 
experiments (table n) are plotted against overall increases in pressure P pk - 1 in 
figure 17, If the data are examined, a Reynolds number effect becomes apparent. For 
each of the three shock generator angles, £/6 0 increases with increasing Reynolds 
number. This condition leads one to seek a consistent Reynolds number effect. 

The initially laminar nature of the boundary layers being considered, the values of 
£/S 0 are divided by \Jr x \ and plotted against P pk - 1 in figure 18. The consistency 
of the results prompts one to seek an exponential variation of the nondimensionalized 
pressure rise lengths with P pk - 1. The resulting equation 

</ 6 o ( p pk - 'f 36 (9) 

,5 

is plotted in figure 18; all the present data fall very close to this curve. 

Equation (9) is empirical and is based only on the present data, which are weak 
incident shock waves interacting with initially ® laminar boundary layers on a sharp flat 
plate for 7 < M 0 < 8, 0.5 x 10 6 < i < 2.0 x 10 6 , and T w *0.5T t< Within this range 
of conditions, equation (9) accounts for boundary -layer thickness, Reynolds number, and 
pressure rise effects. Indeed, the trends of £ with 5 0 , Rj^, and P pk suggested by 
equation (9) are similar to those suggested by other investigators. (See refs. 21 to 25.) 
However, the equation does not account for wall temperature, Mach number, transition, 
and other effects that strongly influence the extent of a laminar pressure rise. (See 
refs. 16 to 25, 33 and 34.) Data from other sources, included in figure 18, do not corre- 
late from the values given by equation (9). The data presented by Watson, Murphy, and 
Rose (ref. 33) were obtained for ^ « 750 000. Their Mach 5.7 data were obtained for 
a ratio of wall temperature to total temperature of approximately 0.41; their Mach 5.8 data 
were obtained for a ratio of wall temperature to total temperature of approximately 0.36. 
The data point from Ball (ref. 21) was obtained for = 330 000 and T w yT t = 0.56. 

The data of Gray and Rhudy (ref. 16) were obtained for adiabatic wall temperatures for 
60 000 < Rx i < 1 000 000. All the data shown in figure 18 are for weak shock interac- 
tions (Pp k < 10j. Data obtained for stronger shocks (such as those presented by Holden 

3 The undisturbed boundary layers are laminar over the flat-plate surface. 
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(ref. 31)) differ greatly from the values indicated by equation (9), depending on wall tem- 
perature, Mach number, and other effects. 

Peak Heating Correlations 

Although many investigations (refs. 5 to 11) have been directed toward obtaining 
correlations of peak heating rates with peak pressures, there is an unacceptably large 
discrepancy in the results for laminar boundary layers. This discrepancy is evidenced 
by a large variation in the values of the exponent n in the simplest, and most frequently 
sought, correlation (refs. 6, 7, 9, and 11): 

Hpk-'Of < 10 > 

Values of n proposed by various investigators are listed in table V . 

TABLE V.- PROPOSED VALUES OF EXPONENT IN CORRELATION 
OF PEAK HEATING WITH PRESSURE (EQ. (10)) 


Investigators 

Suggested n 

Markarian (ref. 6) 

1.29 

Neumann and Burke (ref. 7) 

0.50 

Haslett et al. (ref. 9) 

.70 

Hung and Barnett (ref. 11) 

1.13 


Values of peak heating amplifications obtained during the present experiments (see 
table n) are plotted against measured peak pressure rises in figure 19. Many other data 
points, presented by Hung (ref. 35), are included in the figure. The large data scatter 
apparent in the figure precluded obtaining a reliable, simple correlation in the form of 
equation (10), Hung and Barnett (refs. 11 and 35) have suggested modifying equation (10) 
by using a simple Reynolds number function instead of a proportionality constant. How- 
ever, examination of the present data shown in figure 19 reveals no consistent Reynolds 
number effects. Many investigators have noted large data scatter when attempting sim- 
ple correlations in the form of equation (10) for laminar flows, and have attributed the 
scatter to transition occurring near reattachment of the free shear layer. (See refs. 2, 
8, 9, 11, and 35.) However, transition in the free shear layer depends on the strength of 
the pressure rise as well as on Reynolds number, (See refs. 9 and 26.) Therefore, one 
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would not expect a function of Reynolds number alone to correct satisfactorily the simple 
correlation 4 shown in equation (10). 

Another correlation, suggested by Bushnell and Weinstein (ref. 5), is 

V - fwVslV 0 ' 5 (n) 

^w V e c p \^w Sin< 7 

where p w , V e , and p w pertain to the reattached flow, 5 g ^ is the thickness of the 
shear layer at reattachment, c p is the turning angle of the flow at reattachment, and 
the subscript e refers to conditions at the outer edge of the boundary layer. The data 
presented by Bushnell and Weinstein were measured at the reattachment of separated 
flows on trailing-edge flaps. Keyes and Morris (ref. 10) used the same correlation 
parameters (eq. (11)) for incident shock wave interactions and obtained a larger propor- 
tionality constant than that obtained by Bushnell and Weinstein for the flap reattachment 
data. The correlation presented by Keyes and Morris passes through the present data 
for incident shock interactions, as shown in figure 20. There is considerable scatter 
about the correlation line shown in figure 20. In addition to important transition effects, 
wall temperatures are different for the pressure and heat-transfer data (T w ~0.5Tj. for 
pressure data, T w ~ 0.4T t for heat-transfer data). Although wall temperature effects 
are known to be particularly important for laminar boundary layers (ref. 25), these 
effects could not be investigated by using the present model. 

CONCLUDING REMARKS 

Pressure and heat transfer distributions were obtained on a flat plate in the vicinity 
of incident shock wave interactions with initially laminar boundary layers. The incident 
shock strengths ranged from those insufficient to cause separation to those sufficient to 
cause moderate extents of separated flows. The nominal free- stream Mach number was 
8, and Reynolds numbers, based on plate length, ranged from 0.46 x 10 6 to 2.5 x 10 6 . 

The pressure and heating rate distributions generally follow the established trends 
for nonseparated flows, incipient separation, and laminar free interaction pressure rises. 
The data support previously derived expressions for the minimum pressure rise required 
to cause separation (incipient separation). However, there are no well-developed plateau 
pressure regions in the measured interaction pressure distributions. The upstream 
extent of the pressure rise, divided by the undisturbed boundary- layer thickness, scales 


4 This is contrary to shock wave interactions with turbulent boundary layers; for 

turbulent flows the simple correlation H p k = (P p k) 0 ’ 85 is valid for a wide range of 

flow conditions (ref. 36). ' ' 
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with (Ppb - ^where Pp^ is the peak pressure ratioj and the square root of the 

Reynolds number. 

The heat-transfer measurements were obtained while the wall temperature was 
lower than it was for the pressure measurements. Furthermore, the disturbance caused 
by the incident shock wave usually caused boundary-layer transition prior to the end of 
the interaction region. These factors strongly affect laminar boundary layers. The 
present data do not support previously proposed, simple correlations of peak heating 
with peak pressure. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., November 4, 1974. 
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(a) Model closeup. 

Figure 2.- Photographs of model in Mach 8 variable -density tunnel. 
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(b) Model in injection system chamber beneath test section. 


Figure 2.- Continued. 
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(c) Model in forward position in test section. 


Figure 2.- Concluded 
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Figure 3.- Sketch of model profile. 


















R/m = 1.1 X 10 6 



R/m = 2.0 x 10 6 



R/m = 3.7 X 106 




Distance from leading edge, x, m 

Figure 7.- Interaction pressure ratio distributions for 6 = 1°. 
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R/m = 1.2 x 10 6 



R/m = 1.9 x IQ® 



R/m = 3.7 X 106 



Figure 9.- Schlieren flow photographs for 9 = 3°. 
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Distance from leading edge, x, m 

Figure 10.- Interaction pressure ratio distributions for 9 = 3°. 
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Figure 11.- Heat-transfer coefficient distributions for 0 = 3°. 





R/m = 1.2 x 10 6 



R/m = 2.1 x 106 



R/m = 3.7 x IQ® 



R/m = 6.6 x 10 6 L-74-1167 

Figure 12.- Schlieren flow photographs for 6 = 5°. 
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Distance from leadi ng edge, x, m 

Figure 13.- Interaction pressure ratio distributions for 9 - 5°. 
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Figure 14.- Heat-transfer coefficient distributions for 9 = 5°. 
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Figure 20.- Correlation of peak heating using parameters suggested by 
Bushnell and Weinstein (ref. .5). (See eq. (11).) 



